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ABSTRACT 

Using Hubble Space Telescope (HST) photometry, we characterize the age of the stellar association in the 
vicinity of supernova (SN) 201 Idh and use it to infer the zero-age main sequence mass (Mzams) of the pro- 
genitor star. We find two distinct and significant star formation events with ages of <6 and 11^^ Myrs, and 

+2 



the corresponding Mzams are >29 and 13_f M, 



respectively. These two bursts represent 18;!^% (young) and 
64^\'l% (old) of the total star formation in the last 50 Myrs. Adopting these fractions as probabilities suggests 
that the most probable Mzams is 13;!;j M0. These results are most sensitive to the luminosity function along the 
well-understood main sequence and are less sensitive to uncertain late-stage stellar evolution. Therefore, they 
stand even if the progenitor suffered disruptive post-main-sequence evolution (e.g. eruptive mass loss or bi- 
nary Roche-lobe overflow). Progenitor identification will help to further constrain the appropriate population. 
Even though pre-explosion images show a yellow supergiant (YSG) at the site of the SN, panchromatic SN 
light curves suggest a more compact star as the progenitor. In spite of this, our results suggest an association 
between the YSG and the SN. Not only was the star located at the SN site, but reinforcing an association, 
the star's bolometric luminosity is consistent with the final evolutionary stage of the 17 Myr old star burst. If 
the YSG disappears, then Mzams= 13;l!i ^O' t>ut if it persists, then our results allow the possibility that the 
progenitor was an unseen star of >29 Mq. 

Subject headings: supernovae; general — supernovae: individual (SN 201 Idh) 



1. INTRODUCTION 

Observational measurements of the masses of super- 
nova (SN) progenitors are currently scant and highly un- 
certain. Of the ^20 SNe that have progenitor mass 
constraints, only about half have well-defined masses 
and the rest have upper bounds (lArnett et al 
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1989 



1999 



2002VSma rtt et al.l200aiVan Dvk et al.l2003allbtlS martt et al 
2004; Mau nd et al.ll2005l: iHendrv et al]|2006nLi et al. 2005 
2006, 2007; Gal-Yam et al. 2007; Crockett et al. 2008 
Gal-Yam & Leonard 2009; Smartt et al. 2009; Smartt 2009 
Smith et al.ll201 1 al ibi). In general, stellar evolution theory pre 



diets a clear mapping between the zero-age main sequence 
mass (M/AMs) and the explosion scenario for is o lated massive 
stars (IWoosley et al.l l200l iHeger et al.l 120031; iDessart et all 
120111) . These limited observations suggest t hat the least 
massive stars explode as SN II-P as expected (ISmartt et al.l 
120091) . More massive stars are expected to lose much of their 
mass and explode as H-deficient SNe (lib and Ib/c). How- 
ever, some H-rich SNe (in part icular Iln) h ave been asso- 
ciated with very m assive stars (iGal-Yam & Leonard 20091 
Smith et al. 20 11allibl). F urthermore, theory ( iClaevs et all 
2011; Dessart et al.l 120 111) and the relatively high observed 



rates of H-deficient SNe dSmith et al.l 1201 Ibl) imply that bi- 
nary evolution may figure prominently in producing the H- 
deficient SNe. Even in light of these complications, the max- 
imum Mzams associated with SN II-P is much lower than ex- 
pected (Smartt et al. 2009 ). While tantalizing, these initial re- 
sults are poorly constrained, and even the simple statement 
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that Stars more massive than ^8 M© explode as SNe requires 
more observational constraints. Thus every new supernova 
marks an important opportunity to add a new progenitor mass 
measurement to this poorly-constrained sample. 

Recently, amateurs and the Palomar Transient Factory 
(PTF) collaboration detec ted supernova (SN ) 2011dhinM51, 
the Whirlpool Galaxy (iGriga et al.l 120111; [Silverman et alj 
1201 U lArcavi et al] 1201 Ih . Initial reports classified it as 
II JSilvermanet al] 120111) . but late time spectra reveals 
this explosion to belong to the rare l ib transitional class 
(lArcavi et alJ 120111; iMarion et aP 120111) . indicating that the 
progenitor lost much of its hydrogen envelope due to either 
binary evolution or unknown single-star mass-loss physics. 
The close proxini ity of this SN ((M-ot)o = 29.42, 7.7 Mpc; 
iTonrv et aJ] 1200 IL companion galaxy NGC 5195) and the 
wealth of Hubble Space Telescope (HST) data on M5 1 of- 
fers a rare opportunity to further explore the physics of the 
progenitor 

Using archival HS T photometry, iMaund et alj ( 120111) and 
IVan Dvk et al.l ( 120111) iden tified a yellow s upergi ant (YSG) as 
the progenitor candidate. IVan Dvk et al.l ( 120111) fit the mag- 
nitude and color of the progenitor to evolutionary tracks and 
derive a zero-age-mai n-sequence pr o genito r mass of Mzams= 
18-21 M0. However. JMaund et all (1201 II) argue that an un- 
certain mass-loss history causes the color to be an unreliable 
characteristic of the progenitor Instead, they treat the bolo- 
metric luminosity as an intrinsic property determined by the 
core mass, which is in turn determined by Mzams- Match- 
ing the bolometric luminosity to the very last stages of evo- 
lutionary models, they derive Mzams= 13±3 Mq. To fur- 
ther complicate the situation, recent observations suggest that 
the YSG may not even be the prog enitor. Based upon the 
characteristics of the S N Ught curves. lArcavi et aLJ ( 120111) and 
iSoderberg et al.l (120111) argue that the YSG is not the progen- 
itor and instead suggest a more compact source. 



This discrepancy in mass measurements highlights the fact 
that interpretation of the precursor photometry alone is sen- 
sitive to uncertain mass-loss physics in the final evolutionary 
stages of the models. Furthermore, using precursor imaging 
to measure progenitor masses requires that precursor imag- 
ing exists and that the SN position is known to sub-arcsecond 
precision. The majority of past SNe have neither pre-existing 
HST imaging, nor sufficiently accurate astrometry. A signifi- 
cant fraction of future nearby SNe will lack precursor imaging 
as well. 

In this letter, we address the progenitor mass with an inde- 
pendent approach; we characterize the age of the stellar as- 
sociation surrounding SN 201 Idh and derive the correspond- 
ing MzAMS fo r the progenitor . This t echnique was thoroughly 
described by iGogarten e t al.l (120091) on a SN "impostor" in 
NGC 300 and by Baden es et al.l (I2009h on SN remnants in 
the Magellanic Clouds. In brief, using well-established stel- 
lar population modeling techniques, we can age-date the star 
formation episode that led to the observed SN. The resulting 
age can place strong constraints on the mass of the precursor 
by leveraging the well-understood properties of a large num- 
ber of main-sequence stars. This complementary approach is 
feasible even when there is no precursor imaging, or when the 
SN position is only localized to within several arcseconds. 

In §2, we detail our application of this technique to 
SN 201 Idh, the most distant and only ongoing SN to which it 
has been applied. In §3, we provide the resulting age and mass 
measurements, and in §4, we compare our resulting progeni- 
tor mas s with those ob tained using direct imaging of the pro- 
genitor (iMMOiiLiLbOlJj .\^nDykeLaL20ll)^^_^fe^n^ 
result s to be most consistent with the results of jMaund et alj 
1201 lb and improve upon their uncertainties. 

2. DATA & ANALYSIS 
2.1. Data Acquisition and Photometry 

The region surrounding SN 201 Idh was imaged with HST 
using the Advanced Camera for Surveys (ACS) on January 
20-21, 2005 (Hubble Heritage Team, target name M51-POS5, 
Proposal 10452). For our analysis, we retrieved calibrated 
data from the HST Archive in the deepest two of the four 
available filters, F555W and F814W. Four dithered images of 
exposure time 340 s were available in each filter, providing a 
total exposure time of 1360 s in each band. 

Resolved stellar photometry was performed using the pho- 
tometry pipeli ne developed for the AC S Nearby Galaxy Trea- 
sury program ( iDalcanton et al.ll2009l) . This pipeline uses the 
DOLPHOT stellar photometry package (Dolphin 2000 ) to fit 
the ACS point spread function (PSF) to all of the point sources 
in the images. Fluxes are then converted to Vega magnitudes 
using the standard zero-points and aperture corrections from 
the ACS handbook. To assess photometric errors and com- 
pleteness, at least 10^ fake star tests are performed by insert- 
ing fake stars of known color and magnitude into the data and 
blindly attempting to recover the m with the sam e softw are. 
For details of the quality cuts, see IDalcanton et akl (120091) . 

2.2. Recent Star Formation History Recovery 

For our analysis, we selected stars within 1.4" (^50 pc) of 
the SN location ( Gogarten et al. 2003)- By varying the radius 
from 30 to 100 pc, we verify that 50 pc maximizes the SF in 
the local bursts and minimizes contamination. This resulted 
in 29 stars being selected, including 16 upper main sequence 
stars. A color-magnitude diagram (CMD) of these stars is 
displayed in the left panel of Fig. [T] 



We fit the stellar evolution models of iGirardi et alj ( 120021 
120101) to our data using the software package MATCH, which 
calculates a maximum likelihood fit of a linear combination of 
model CMDs to the observed CMD (Dolphin 2002J. We veri- 
fied that our CMD was consistent with a populati on located at 
the pu blished distance of M51 ((m-M)n=29.42. lTonry et alJ 
1200 II) by allowing MATCH to fit the distance modulus. The 
best-fit distance modulus is (ot-M)o=29.45 and is consistent 
with the published distance. MATCH also fits for a distribu- 
tion of extinctions; we fix the width of the distribution at 0.5, 
but allow the minimum to vary from Av=0.0 to 0.5. We con- 
strained the metallicity to increase for younger ages. We only 
fit data brighter than our 50% completeness limit as measured 
by our artificial star tests. These limits were m <26.95 in both 
filters. 

Having calculated a best-fit age distribution, we quantify 
the random and systematic uncertainties. The random error 
is due to sampling statistics and is largely determined by the 
number of stars in the sample. The systematic error is due 
to potential offsets between model and observed CMDs. For 
example, if the model is systematically slightly bluer than the 
data (potentially due to systematic offsets in the model trans- 
formations or model deficiencies), then the resulting age dis- 
tribution will be affected. We estimate the uncertainty from 
both of these sources through a series of Monte Carlo (MC) 
realizations of the data. Each realization draws the same num- 
ber of stars as our data from the best-fit age distribution. 
These random draws are re-fit with the models after shift- 
ing the models by 0.4 magnitudes in bolometric luminosity 
and 7% in effective temperature, which represent distance and 
systematic model uncertainties. Thus, the set of resulting fits 
provide a distribution of measurements from which we calcu- 
late our uncertainty. Finally, to assess the dependence of the 
result on individual outliers, we performed jackknifing tests 
in which the brightest or bluest star in the field (such as the 
SN progenitor candidate) was removed from the data. 

Our goal in this analysis is to identify the burst of star for- 
mation in which the progenitor star was most likely created. 
We therefore focus on star formation history (SFH) of the past 
50 Myr To identify recent star formation bursts, we examine 
the cumulative star formation as a fraction of the total stellar 
mass formed in the past 50 Myr We then calculate the median 
age of the recently formed stars as well as the width of their 
age distribution. This method has the advantage of only being 
sensitive to relative star formation rates, not absolute values. 
Furthermore, the cumulative distribution is less affected by 
the covariances between neighboring age bins in the model 
fits. 

This methodology leads to two sources of uncertainty in 
age. One is the uncertainty in the median age, which the MC 
tests address, and the second is the width of the distribution of 
ages. To account for the uncertainty due to the intrinsic age, 
we assign probabilities to each age based on the mass fraction 
of stars with that age. Uncertainties on these probabilities 
then account for the measurement errors. Our quoted mass 
values correspond to the most massive star for the age limits, 
which are calculated directly from the same models as used 
for fitting (IGirardi et al.ll20IOl) . 

3. RESULTS 

In this section, we report the best-fit metallicity, average 
extinction, and SFH of the region surrounding the location of 
SN 201 Idh. In Fig. [2] we show the primary result, the age dis- 
tribution of young populations surrounding SN 201 Idh. The 



best-fit metallicity for the last 50 Myr is solar, which com- 
pares favorably with the range of metalUc ities derived from 
HII regions in M51 ( iBresolin et al.ll2004l) . For the average 
extinction, we find Ay = 0.425. The value for galactic extinc- 
tion at this pointing is Ay=0.12 (iSchlegel et al.lll998l) : there- 
fore this region likely contains additional extinction from dust 
within M51. To estimate the uncertainty due to extinction, 
we force A,,=0.12 in MATCH and compare these results to 
the self-consistent fit. The resulting differences in ages and 
masses are included in the uncertainties that we report in this 
letter. 

The left panel of Fig. [T] shows the magnitude (F555W) 
and color (F555W-F814W) of the brightest stars surround- 
ing the SN, which we use to calculate the SFH in the vicin- 
ity of SN 2011dh. The green star shows th e magnitude 
and color of the progenito r candidate, a YSG (iMaund et al.l 
1201 U IVan Dyk et al.l 120 llh . For comparison, we plot three 
isochrones corresponding to the three star formation (SF) 
events in Fig. |2] In addition, we show the corresponding age 
and fraction of total SF in the last 50 Myrs. The right panel 
of Fig. [U shows a model CMD produced by MATCH for the 
best-fit SFH, with the same isochrones over-plotted. Note that 
the model generically reproduces the number density of stars 
on the MS and red giant regions of the CMD. 

To be clear, we did not fit the main sequence turn-off 
(MTO) to arrive at these isochrones. Because of the small 
number of main sequence stars at the bright end, MTO fitting 
is susceptible to large Poisson errors. Rather, as we explain in 
§ |2l we model the entire CMD and constrain the SFH that is 
required to reproduce the observed CMD. This way the num- 
ber density of stars at each color and magnitude, particularly 
along the entire MS, is involved in constraining the age of the 
star burst, not just the location of the end of the MS. 

The best fit SFH is shown in Fig. |2] The left panel shows 
the relative SFH (black solid line) as a function of age, and 
the right panel shows the cumulative SF (thick solid purple 
line) since 50 Myr ago. The orange thin lines show the 68% 
confidence interval for the cumulative distribution. In both 
plots, three bursts of SF are apparent: 4, 8, and 17 Myr The 
17 Myr burst accounts for most (64;!;[4%) of the recent SF in 
the last 50 Myr 

Based on these bursts, we report two estimates for the age 
and MzAMS- The age is either <6 or ll'!il Myr, corresponding 
to MzAMS of either >29 or 13;I;i ^0' respectively. Since most 
of the SF in the last 50 Myr is associated with the ~17 Myr 
burst, we find the most probable Mzams of the progenitor to 
be 13;!;j M0. We rule out the ~8 Myr SF event for three inter- 
related reasons. For one, this event has uncertainties con- 
sistent with zero star formation. Secondly, jackknifing tests 
show that this event is strongly dependent on the presence 
of the YSG, suggesting that no other stars in the sample are 
consistent with this age. Thirdly, if the YSG is the progeni- 
tor, then its color implies that it experienced a great deal of 
mass-loss that is not included in the models used in MATCH. 
Hence, any fits that rely on the magnitude and color of this star 
alone are in error. See §|4]for further discussion on these last 
two points. In the cumulative plot, the age and mass estimate 
associated with the younger event is highlighted by the blue 
shaded regions, while the older and most probable estimate is 
highlighted by the red shaded region. 



4. DISCUSSION & CONCLUSIONS 



Using HST photometry, we characterize the age of the stars 
in the vicinity of SN 201 Idh, and from this age, we infer the 
A/zAMS of the progenitor star The recent SFH shows two 
SF events with ages of <6 and ll'^i^ Myrs, corresponding to 
A^ZAMS for the progenitor of >29 and 13;I;i ^O' respectively. 
Given that the older burst at 17;!:4 Myr represents 64;!: [4% of 
the total SF in the last 50 Myrs, the most probable Mzams for 
the progenitor is 13;!;[ M0. Because MATCH leverages the 
entire CMD, including the MS, these MzAMsestimates are rel- 
atively insensitive to the individual peculiarities of post-MS 
or binary evolution. 

Stellar population analysis shows that the colors of the YSG 
are consistent with the surrounding stars only if it experi- 
enced peculiar evolution. Testing the robustness of the three 
SF events, we perform several jackknifing tests, removing a 
bright star and finding new best-fit SFHs. The <6 and 17 Myr 
old populations are robust to these tests, but upon removing 
the YSG, the 8 Myr old burst vanishes (<3.5% of the recent 
SF). Hence, the 8 Myr feature is solely dependent upon the 
color and magnitude of the most poorly modeled star in the 
field, the YSG - an unlikely scenario. More likely, the YSG 
was born in one of the other events and experienced disruptive 
mass loss, giving it peculiar colors. In fact, its bolometric lu- 
minosity is consistent with the final evolutionary stage of the 
17 Myr old population. 

Regardl ess of an association with the SN, the lum inosity of 
the YSG (IMaund et all 120111: 1 Van Dyk et all 1201 Ih suggests 
that it is associated with the 17 Myr population, in which case, 
we derive a Mzams for the YSG of O;!!; Mq. The previous 
Mzams estimates were obtained by comparing photometry of 
the YSG with stellar evolution models. Of the two previous 
attempts to estimate Mzams one at 13 ±3 Mq (|Maund et al.l 
120111) and one at 18-21 Mq (IVan Dyk et al.ll201 ll) . the lower 
mass estimate is most consistent with our results. This con- 
sistency leads to a couple of conclusions. 

For one, this consistency validates using nearby stellar ages 
to estimate the progenitor mass when an archival HST im- 
age of the progenitor is not available. Given the distance to 
M51, ~7.7 Mpc, the fact that this technique is able to further 
constrain the Mzams estimate is remarkable. Directly model- 
ing the progenitor in precursor imaging only requires accurate 
photometry for the most luminous stars and has become a rou- 
tine exercise up to 20 Mpc ([Smartt 2003). Hence, for direct 
imaging of the precursor, SN 201 Idh is quite close. Our tech- 
nique, on the other hand, relies on having a sufficient number 
of main sequence stars to characterize the age of the associ- 
ated star burst. With a distance of ^7.7 Mpc, the magnitude 
limit is quite high, resulting in only 16 detectable upper main 
sequence stars within 50 pc of the SN. Even with this small 
number of upper main sequence stars, we were able to im- 
prove upon the direct imaging mass constraint for the YSG. 

The second conclusion we draw from the consistency is that 
fitting the bolometric lumino sity of the progeni tor to the last 
stages of evolutionary tracks (IMaund et al.l201 11) is a more ro- 
bust method to estimate Mza ms than fitting both th e luminos- 
ity and effective temperature dVan Dyk et al.ll20Tlb . The liter- 
ature is full of theoretical arguments suggesting that the bolo- 
metric luminosity during the last evolutionary stages is a more 
intrinsic property of Mzams (lArnett et al.ll989l:IWoosley et al.l 
1^02; Smartt 2009), but as far as we know, there have been no 
independent observations to support this. Because age-dating 
the nearby stars leverages information from the MS phase, our 
results provide a complimentary estimate of Mzams and inde- 



pendent support for using the bolometric luminosity only in 
direct imaging techniques. 

While the most probable Mzams of the SN progenitor is 
13!!;[ M0, there is a small but not insignificant probability that 
the progenitor mass is >29 M©. Ruling out one or the other 
as the progenitor requires further information. If the YSG per- 
sists i n post-SN irnaging , then the progenitor wa s indeed com- 
pact (lArcavi et al.ll201 U [Soderberg et"aI1l201 lb . and neither 
mass estimate is ruled out. However, if the YSG disappears, 
then our progenitor mass esti mate of 13;!;? M w would validate 
and further constrain the Maund et al.l (1201 lb estimate. Given 
the YSG's positional coincidence with the SN and the fact 
that its bolometric luminosity is consistent with the final evo- 
lutionary stage of the 17 Myr old population, it would be odd 
if the YSG is not associated with the SN. 
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Figure 1. Observed (left panel) and modeled (right panel) c olor-magnitude diagrams (CMD). In the left panel, the brightest star (green star symbol) is the 
progenitor candidate jMaund et a l.ll201 llJVan D vk et alj|201 II ). The three lines correspond to Padova isochrones that are closest in age to the three SF bursts in 
Fig!l2]and each is labeled with the age and percent of total SF in the last 50 Myr. Note that we do not fit these isochrones to determine the age. Rather, we model 
the CMD and find the SFH that best represents the data. The right panel shows a representative CMD model from the MC simulations. 
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Figure 2. Modeled star formation history (SFH) in the vicinity of SN 201 Idh. The left panel shows the relative star formation rate (solid black line). We use the 
cumulative SF since 50 Myr ago (right panel) to characterize the age. We find two significant SF events at <6 Myr and 17^;^ Myr (See §|4]for a discussion of 
how we rule out the 8 Myr old burst). The corresponding Mzams for the progenitor are >29 and 13;';^ Mq. Most of the recent SF occurred in the 17 Myr burst. 
Therefore, our most probable Mzams estimate for the progenitor is 1 3+^ Mq . 



